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Abstract Using a recent magnetohydrodynamic simulation of magnetotail dynamics, we further 
investigate the buildup and evolution of the substorm current wedge (SCW), resulting from flow bursts 
generated by near-tail reconnection. Each flow burst generates an individual current wedge, which includes 
the reduction of cross-tail current and the diversion to region 1 (Rl)-type field-aligned currents (earthward 
on the dawn and tailward on the duskside), connecting the tail with the ionosphere. Multiple flow bursts 
generate initially multiple SCW patterns, which at later times combine to a wider single SCW pattern. The 
standard SCW model is modified by the addition of several current loops, related to particular magnetic field 
changes: the increase of B z in a local equatorial region (dipolarization), the decrease of \B X \ away from the 
equator (current disruption), and increases in \B y \ resulting from azimuthally deflected flows. The associated 
loop currents are found to be of similar magnitude, 0.1 -0.3 MA. The combined effect requires the addition of 
region 2 (R2)-type currents closing in the near tail through dawnward currents but also connecting radially 
with the R1 currents. The current closure at the inner boundary, taken as a crude proxy of an idealized 
ionosphere, demonstrates westward currents as postulated in the original SCW picture as well as 
North-South currents connecting R1 - and R2-type currents, which were larger than the westward currents 
by a factor of almost 2. However, this result should be applied with caution to the ionosphere because of our 
neglect of finite resistance and Hall effects. 


1. Introduction 

The substorm current wedge (SCW) [ McPherron et al., 1973] constitutes the main paradigm of 
magnetotail-ionosphere interaction associated not only with substorms but also with earthward flow burst 
events that do not lead to full substorms [e.g., Henderson et al., 1 998; Sergeev et at, 1 999; Lyons et al., 1 999; 
Nakamura etol, 2001a, 2001b]. Its major component consists of a reduction of the cross-tail current in a cen- 
tral tail section, associated with a collapse and dipolarization of magnetic field, combined with a deviation 
of current to the ionosphere, where it closes through the westward auroral electrojet. This simple current 
loop involves field-aligned currents connecting magnetotail and ionosphere with the sense of "region 1" 
(R1) [Hjima and Potemra , 1976], that is, earthward on the dawnside and tailward on the duskside. 

In a recent paper [ Birn and Hesse , 201 3] we used magnetohydrodynamic (MHD) simulations to demonstrate 
the buildup of the field-aligned current system by flow bursts or "flow channels," shearing, or twisting the 
magnetic field. These simulations also showed the relevance of close penetration of the flows in the buildup, 
the incorporation of an additional, region 2 (R2) type, current system, and the possibility of multiple flow 
bursts leading to an accumulation of dipolarized flux tubes. 

Using the same MHD simulation of near-tail reconnection [Birn et al., 201 1 ], we here illustrate further proper- 
ties of the current wedge picture, focusing in particular on the different current systems and their evolution. 
The description of the SCW, or any current system, by some number of current loops appears somewhat 
arbitrary. Here we demonstrate that a natural way of identifying major current systems can be associated 
with the main changes in the magnetic field components. Major properties of the MHD simulation are 
described by Birn etal. [2011] and Birn and Hesse [2013] and will not be repeated here. Here we just note that 
the simulation was based on dimensionless quantities with suitable units given by 

L n = 10,000km « 1.5 R E , B n = 20nT, v„ = 1,000km/s (1) 

This leads to derived units t n = L n /v n = 10 s, p n = B 2 n /pi 0 = 0.32nP,y n = B n /(pi 0 L n ) = 1.6nA/m 2 , and 
l n = B n L n /fi 0 = 0.16 MA. The simulation box spans the region -7.5 R E > x > -97.5 R E , \z\ < 15 R E , 
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Figure 1 . Entropy function S = f p } / y ds/B (color) and flow vectors in thex,y plane at three different times, for a simu- 
lation of near-tail reconnection, modified after Figure 6 of Birn and Hesse [201 3]. Solid black lines represent contours of 
constant B z , shown at intervals of 0.5 (10 nT); the one farthest to the right represents the B z = 0 line. 


|y| < 60 R e . (The extent in y was unusually wide for comparison with a simulation with strong cross-tail mag- 
netic field, which is not relevant here.) The initial state consisted of a tail field with a superposed dipole with 
a center located at x = 0 outside the simulation box. The evolution included a relaxation phase (0 < t < 30, 
corresponding to 300 s), during which the system relaxed into full equilibrium [Hesse and Birn, 1993], fol- 
lowed by a driven phase (30 < t < 61 ), during which an external inflow of magnetic flux was applied to the 
top and bottom boundaries. This leads to current intensification and the formation of a thin embedded cur- 
rent sheet. At t = 61 finite resistivity was imposed, concentrated in the region of enhanced current density, 
which caused the onset of reconnection and the formation of a neutral line (B z = 0) at t « 90. 

A main property of interest for the present paper is the fact that the simulation consisted not only of an 
initial flow burst near midnight, starting around t = 125, but also of additional flow bursts off midnight 
delayed by a few minutes. These flow bursts and their association with entropy-depleted magnetic flux 
tubes are illustrated in Figure 1, modified and extended after Birn and Hesse [201 3]. Here the color-coded 
function 5 denotes an entropy measure defined by 

5 = J p 1/7 ds/8 with y = 5/3 (2) 

integrated along field lines from one boundary crossing point to another [Birn etal . , 2009; Wolf eta!., 2009]. 
For closed field lines of interest here the crossing points are at the near-Earth boundary. (For plasmoid field 
lines, identified by several crossings of the z = 0 plane, the integral was taken over one loop, from one 
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equatorial crossing to another.) Solid 
black lines represent contours of con- 
stant B z , shown at intervals of 1 0 nT; 
the one farthest to the right represents 
the£ z = 0 line. The closed contours at 
the front of each flow burst indicate the 
dipolarization fronts. Black arrows show 
velocity vectors with a speed above 
100 km/s. 

At t = 1 30, Figure 1 shows a single flow 
burst near midnight, associated with a 
region of reduced entropy originating 
from the reconnection site. The velocity 
vectors demonstrate not only the earth- 
ward flow bursts but also flow vortices 
surrounding the dipolarization front. At 
t = 1 60 the central flow burst has largely 
subsided, leaving a region of enhanced 
B z and reduced 5 near the dipolar region, 
while additional flow bursts neary = ±7 

can be seen, each associated with their individual low-entropy, high-£ z regions, and vorticity patterns. At 
t = 200 these flow bursts have been deflected and slowed, while additional flows within |y| < 3 add to the 
dipolarization near midnight. Finally, at t = 250 the near-Earth flows have largely subsided and left a wide 
accumulation of low-entropy, high-£ z regions. 



Figure 2. Evolution of the total current of region 1 (positive, red curve) 
and region 2 signature (negative, green curve), integrated fory < 0,z > 
0 within the closed field line region. The current unit corresponds to 
0.16 MA. 
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Figure 3. Color-coded field-aligned current component j^ x at the inner 
boundary for various times as indicated. The dashed lines show the 
open-closed field line boundaries. 


Section 2 shows the evolution of the 
current systems generated by the flow 
bursts, while in section 3 we investigate 
more closely the various current sys- 
tems associated with the magnetic field 
changes driven by an individual flow 
burst, which results in a more complete 
SCW picture. Section 4 addresses the cur- 
rent closure at the inner boundary, taken 
as a crude proxy for the ionosphere, 
and section 5 provides a summary 
and conclusions. 


2. Evolution of Current Systems 

As demonstrated by Figure 1, the ini- 
tial flow burst near midnight is followed 
by additional flow bursts off midnight, 
leading to a somewhat turbulent accu- 
mulation of magnetic flux nearer to 
Earth. Each of the flow channels gen- 
erates current systems involving also 
field-aligned currents that connect the 
tail with the ionosphere. First, we pro- 
vide an overview of the evolution of 
the total region 1- and region 2-type 
field-aligned currents at the inner simu- 
lation boundary, shown in Figure 2. The 
curve for the region 1-type current (red 
curve) differs slightly from Figure 7b in 
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t = 250 



Figure 4. Color-coded divergence to field-aligned currents V • Iy, with Iy = / jydz, together with changes of the current 
density vectors at z = 0, Aj = j - j 0 . Solid lines are contours of constant B z at intervals of 0.5 (1 0 nT). 


Birn and Hesse [2013], because here we included only currents within the closed field line region, whereas 
Birn and Hesse [201 3] also included currents within the open (lobe) regions. Figure 2 also includes the driven 

phase (30 < t < 61, dashed lines), which 
already causes moderate increases in both, 
region 1 and region 2, systems. The integrated 
effects of the major flow channels at t > 125 
and t > 150 are clearly visible, resulting in 
increases of both region 1- and region 2-type 
currents. The increases in region 1-type cur- 
rents are about twice as large as those of the 
region 2 currents. We note that these currents 
most likely do not represent the total current 
closing within the ionosphere, because they 
may be comprised of individual "wedgelets" 
[Rostoker, 1 998], which may close separately. 

Details of the field-aligned currents at the inner 
boundary are demonstrated by Figure 3, which 
shows the color-coded field-aligned current 
component y‘y x for various times, demonstrat- 
ing the effects of the initial flow burst (f = 140) 
and the two subsequent flow bursts (f = 1 60), 
each causing a wedgelet of region 1 -type with 
a minor region 2-type system at lower latitude. 
At later times the effects of the flow bursts 
decrease within the expanding bulge, while the 
Figure 5. Schematic of major current systems contributing to major contributions from the edges remain, 
the SCW, modified from Plate 4 of Birn et al. [1 999] by the addi- 

tion of loops 4 and 5; a preliminary version of this figure without The evolution of current diversion (color coded) 
loop 5 is also shown by Kepko et al. [2014], and the changes in the currents (vectors) in the 
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20 
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Figure 6. (a) Divergence to field-aligned currents, integrated 
over z (color), flow vectors, and magnetic field ( B z ) contours in 
the equatorial plane at t = 132; (b) equatorial current density 
vectors Aj (arrows) associated with dipolarization (A B z , color), 
indicating loops 4 and 5 in Figure 5; blue and yellow contours 
show the regions of main current diversion from Figure 6a. 


equatorial plane are demonstrated by Figure 4. 
Here A j = j - j 0 are the changes of the current 
density vectors at z = 0, where j 0 represents 
current density vectors at t = 61 prior to the 
onset of reconnection. While the effects of the 
initial flow burst remain near y = 0, the two 
additional flow bursts neary = ±7 (Figure 1) 
each sets up an additional current loop of 
type 4 (Figure 5) with the corresponding cur- 
rent diversion predominantly to region 1. 
However, at late times the current eddies within 
the inner region disintegrate, and the system 
appears dominated by a single wide loop 
of type 4. 

3. Modified SCW Loop Picture 

In this section we look more closely into 
the current systems set up by an individ- 
ual flow burst. We start with a schematic of 
the major current systems that we found 
to contribute to the effects of a single flow 
burst. These current systems are illustrated 
in Figure 5, modified from Birn etal. [1999, 

Plate 4] by the addition of loops 4 and 5. (An 
earlier version of this figure without loop 5 
is presented by Kepko etal [2014]; Figure 3.) 
The presence of loop 4 has also been demon- 
strated earlier in Plate 1 of Birn etal [2000]. 


AB X 
■“ 1 .0 


Loop 4 is directly related to the increase of B z 
near the equatorial plane, that is, the dipolar- 
ization, as demonstrated by Figure 6b, showing 

the changes in current density, A j = j - j 0 (where j 0 again represents current density vectors at f = 61 
prior to the onset of reconnection), superposed on the color-coded magnitude of A B z = B z - B z0 . The 
top part shows the color-coded regions of current diversion V • Iy, obtained by integrating V • jy in z, 

together with contours of B z and velocity 
vectors. The regions of current diversion are 
also indicated in Figure 6b by heavy blue and 
yellow contours. 

Figure 6a shows not only the dominant diver- 
sion to region 1-type currents, resulting from 
the vortex flows indicated by the velocity 
vectors, but also a diversion to region 2-type 
currents earthward of it, pointed out already 
by Birn and Hesse [201 3, Figure 1 5]. This sys- 
tem is set up by the flow shear earthward of the 
azimuthally deflected flow. Figure 6b indicates 
that this diversion is spatially related to current 
loop 5 in Figure 5 with opposite sense to loop 
4, surrounding a region of slightly decreased 
magnitude of B z . This region coincides 
with a region of enhanced pressure [Birn and 
Hesse , 201 3, Figure 9] just earthward of the 
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Figure 7. Changes in B x (color) at t = 132 for x = -12, together 
with current density vectors (Aj), related to loop 1 and 3 cur- 
rents in Figure 5. The colored contours indicate the regions of 
enhanced parallel currents. Dashed lines show the open-closed 
field line boundaries. 
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dipolarization front. The current diver- 
sion is related to azimuthal pressure 
gradients, as governed by Vasyliunas' 
equation [ Vasyliunas , 1970] 


j\\ = — (B/8) • Vp X VI/ (3) 

where V = J ds/B is the differential 
flux tube volume, integrated along a 
field line. 


Figure 8. (a) Color-coded B y and (b) v y at t = 1 32 at y = -1 .5, together 
with current density vectors (Aj). Heavy blue and orange contours indi- 
cate regions of enhanced tailward and earthward field-aligned currents, 
respectively. Dashed lines show the separatrices (open-closed field 
line boundaries). 


Another effect of the dipolarization is the 
decrease of |8 X | at higher latitudes off 
the equatorial plane, demonstrated by 
Figure 7, which is closely related to the 
reduction of the total current (current 
disruption). This reduction is associ- 
ated with loops 1 and 3. While loop 1 
represents the main part of the stan- 
dard SCW picture, loop 3 indicates that 
part of this current may close within 
the magnetotail, which is presum- 
ably the dominant type in Mercury's 
magnetosphere, where the lack of an 
ionosphere prevents the possibility of 
ionospheric current closure [Kepko 
et al., 2014]. 


The additional current systems shown 
as loops 2 are related to the braking 
and diversion of flow around the Earth, 
which results in an increase of \B y \ , 

illustrated in Figure 8a for the dawnside by red (z > 0) and blue (z < 0) colors, surrounded by earth- 
ward current at higher latitude and tailward current at lower latitude. The connection between these 
currents occurs away from the equator, while near the equatorial plane the earthward currents of loop 
4 dominate. 


AJ V 




Figure 8a also shows regions of opposite B y at higher latitudes (-13 < x < —10, |y| « 3) just inside the 
separatrices (dashed lines). As demonstrated by Figure 8b, these distortions of B y result from flows toward 
midnight (red-orange regions), which are part of the vortices generated on the outside of an earthward 
flow burst (Figure 1 ). The distortions of B y lead to field-aligned currents, opposite to those of region 1 , but at 

higher latitude. In the simulation, 
however, these currents are mainly 
closed within the box and do not lead 
to notable field-aligned currents at 
the inner boundary. For simplicity, 
these current loops are not included in 
Figure 5. 

All loops, except for loop 2, close through 
the midnight meridional plane, y = 0. 
This allows us to estimate their relative 
magnitude. Figure 9 shows the 
color-coded change in the cross-tail cur- 
rent density at t = 132, red showing 

Figure 9. Color-coded current density Aj y at f = 1 32 at y = 0. Red and duskward and blue dawnward cur- 
blue colors show duskward and dawnward currents, respectively. The rents. The currents are by no means 

dashed lines emphasize the boundaries between the two directions. "wirelike." They form layers or are even 
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Table 1. Magnitude of the Cross-Tail Currents in Regions l-V in Figure 9 and the Resulting Currents in Loops 
1-5, Schematically Shown in Figure 5, at t = 132 

Loop /, /„ In, l/ V l V 1 2 3 4 5 R1 R2 

A/(MA) 0.31 0.51 0.69 0.59 0.47 0.37 0.21 0.12 0.20 0.31 0.58 0.21 


volume-filling regions in space. The blue region earthward of x > -9 (region I) corresponds to the dawn- 
ward portion of loop 5. The adjacent red region earthward ofx > -10 (region II) consists of the duskward 
parts of loops 5 and 4. The blue region further out, extending to x « -18atz = 0 (region III), consists of 
the dawnward parts of loops 1 , 3, and 4. This region is surrounded by thin red layers connected to the recon- 
nection site (region IV), which include the duskward closure of loop 3. These layers are surrounded at larger 
|z| by dawnward currents (blue, region V), which presumably close with parts of the underlying duskward 
currents. (This loop is not included in Figure 5.) 

We integrated the eastward and westward currents in Figure 9 for z > 0 out tox = -18 to provide estimates 
of the relative magnitudes of the different loops within one hemisphere. The results are shown in Table 1 . 
The loop 2 current was assumed to be ~ 60% of the loop 1 current, based on our finding that the change in 
total region 2 (R2)-type current at this time is about 37% of the change in region 1 (Rl)-type current (com- 
prising both loop 1 and loop 2 currents), as indicated by Figure 2. This ratio agrees well with observational 
results [ Sergeev et al, 2014]. 



Figure 10. Magnetic distortion from (a) dipolarization and (b) 
azimuthal flow. 


4. Ionospheric Current Closure 

On the short substorm expansion time scale 
we might, to the lowest order, consider mag- 
netic field lines as tied in the ionosphere, that is, 
treating the ionosphere as a solid, ideally con- 
ducting wall. This is equivalent to the boundary 
conditions imposed at the inner boundary 
of our simulation box. The closure of the cur- 
rents at this boundary then might provide a 
crude proxy for the ionospheric closure. The 
boundary currents can be considered as sur- 
face currents associated with the perturbations 
of B z and B y across the conducting surface with 

Jy = f jydX = ~(B Z+ - B z _ ) (4) 

where B z _ and B z+ define the values of B z just 
inside and outside the inner boundary, respec- 
tively. If we assume that the value outside 
remains unchanged, we can identify B z+ with 
the value at the initial time, t = 61 , and obtain 

J y = AB z = B z -B z \ t ^ (5) 

Similarly, we find 

J z = f j z dx = (B y+ - B y _) = -AS y (6) 

Thus, positive A B z represents westward current. 
(Note that this corresponds to perturbations 
above the ionosphere.) Positive (negative) 

A B y represents equatorward (poleward) cur- 
rent on the dawnside (duskside) Northern 
Hemisphere with opposite directions in the 
Southern Hemisphere. 
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Figure 1 1 . Magnetic perturbations, A B z and A B y (color), at the inner boundary for various times as indicated. Arrows 
indicate the surface current vectors associated with these perturbations. The colored contours indicate the regions of 
strong field-aligned currents j^ x shown in Figure 3. The dashed lines again show the open-closed field line boundaries. 


The basic relationship between the magnetic field perturbations and the currents is illustrated in Figure 1 0. 
Figure 10a shows a dipolarized field line (red) between two still-stretched field lines (blue). The shear 
between the red and the blue lines (rotation of B z /B x as function of y) is the signature of the (region 1 -type) 
field-aligned currents and the increase of B z at the boundary indicates westward current at the foot points 
of the affected field lines. Figure 1 0b shows an azimuthally bent field line (green) between two field lines 
(blue) more closely aligned with the x direction. The shear between the green and the blue lines (rotation of 
B y /B x as function of z) is the signature of the field-aligned currents of region 1 signature above and region 
2 signature below the distorted green field line; the increase of \B y \ at the boundary indicates equatorward 
current at the foot points of the affected field lines. 

The current closure in the simulation is demonstrated clearly by Figure 1 1, which shows the color-coded 
perturbations of the magnetic field components A B z and A B , perpendicular to the magnetic field, and 
the associated surface current vectors (arrows) at the inner boundary for various times. Regions of strong 
field-aligned currents j^ x , as shown in Figure 3, are indicated by colored contours and the open-closed 
boundaries by dashed lines. The perturbations A B z (Figure 1 1, left) represent current closure through west- 
ward currents, which is part of the traditional SCW picture, whereas the perturbations A B y (Figure 1 1, right) 
show the North-South closure between the region 1 - and region 2-type currents. The integrated currents at 
t = 140 correspond to ~2.6x 10 5 A (westward) and ~4.6 x10 5 A (equatorward on the dawnside). It is not 
clear, however, to what extent this result can be transferred to ionospheric closure, because of our crude 
treatment of the "ionosphere" as a plane surface, neglecting finite resistance and Hall effects, which are 
expected to contribute to an enhanced closure through east-west currents through the "Cowling channel" 
[e.g.. Cowling, 1 932; Baumjohann, 1 983; Yoshikawo et ol., 2013]. 

At t = 160 each of the wedgelets is associated with both types of current closures; there are only minor 
eastward currents between different wedgelets, whereas at later times (f = 200) the current distribution 
approaches a single wider wedge. The westward currents at this time, however, are not very uniform and 
spread over different latitudes. Showing the perpendicular components of A B z eliminates compressional 
effects near z = 0, which are not expected to contribute to ionospheric current closure. 
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5. Summary and Discussion 

Using a three-dimensional MHD simulation of magnetotail reconnection and dipolarization [ Birn et al., 
2011], we have further investigated properties of current diversion associated with the substorm current 
wedge (SCW), extending the investigation of Birn and Hesse [201 3]. Particular emphasis was on the current 
systems contributing to the SCW during various stages. In this context it is useful to distinguish the buildup 
of the current systems and their persistence. 

In general, the buildup of the SCW can be attributed to the consequences of one or more earthward 
flow bursts, generated by reconnection and apparently modified by cross-tail structure associated with 
ballooning/interchange-type instability [Birn et al. , 201 1 ]. The evolution starts with a single flow burst near 
midnight, followed by two flow bursts off midnight and subsequent more turbulent flows. The flow chan- 
nels can be identified as low-entropy flux tubes, associated with strong enhancements of B z (dipolarization 
fronts). The low-entropy, high-£ z flux tubes accumulate near Earth when the flows get stopped and diverted 
around Earth. Each of the flow channels distorts the magnetic field and builds up (drives) a current system 
akin to the standard SCW (a wedgelet) [ Rostoker , 1998], but when the dipolarized fields accumulate near 
Earth the combined system eventually becomes dominated by a single wide wedge. 

The distortions of the magnetic field and the associated currents persists even after the flows subside. In 
that case the current systems and particularly the field-aligned currents are commonly inferred from force 
balance [e.g., Wolf et al., 2009], involving pressure gradients and gradients of the magnetic flux tube volume. 
Here we focused instead on the magnetic field changes, which are more easily observable than the cur- 
rents and can be associated with simple current loops, supporting and complementing the standard SCW 
picture. The major current systems contributing to the total SCW as well as to a single wedgelet are illus- 
trated in Figure 5, modified from Birn etal. [1999]. They can be attributed to localized changes in B x , B yl and 
B z , respectively. 

One main effect of a flow burst is the dipolarization, which is part of a propagating low-entropy flow chan- 
nel but also part of the accumulation of deposited low-entropy flux tubes. The enhanced B z at or near the 
equatorial plane is associated with current loop 4, reducing the cross-tail current farther out but increasing it 
closer in with radial connections at the edges of the dipolarized region. A minor modification is current loop 
5, which is related to a local reduction of B z , associated with a pressure increase earthward of a dipolarization 
front, with azimuthal gradients that correspond to a diversion to region 2-type currents. 

A second effect of the dipolarization is the decrease of \B X \ at higher latitudes, which is closely related to 
the reduction of the total current in that region (current disruption). This reduction is associated with loops 

1 and 3. These loops contribute to the region 1-type field-aligned currents, which are part of the standard 
SCW. In Mercury's magnetosphere, where the lack of an ionosphere prevents the possibility of ionospheric 
current closure [ Kepko et al., 2014], they presumably close within the magnetotail (loop 3). 

The third effect is the increase of \B y \ resulting from the field deformation by azimuthal flows after braking 
and diversion. This results in the current loops 2 contributing to both region 1 on the outside and region 

2 currents on the inside. The flow braking is also associated with a pressure enhancement, which provides 
an alternative method of inferring the region 2-type currents via equation (3) from the azimuthal pressure 
gradients. This result is consistent with Time History of Events and Macroscale Interactions during Substorm 
(THEMIS) observations [ Dubyagin et al., 2010]. 

The currents in loop 4 in Figure 5 appear to be unrelated to the wedge currents. However, a different view 
arises when the currents of Figure 5 are combined, disregarding loops 3 and 5, which are confined to the 
tail. In that case we obtain the two basic loops of Figure 12, equivalent to Figure 1 b of Sergeev et al. [2014]. 
Part of the loop 4 currents now contribute to the region 1 currents, while the earthward portion becomes 
part of the region 2 loop. We note, however, that the region 1 and 2 loops are not separated but may be con- 
nected through radial tail currents (Figure 1 2, dashed lines). These currents were found to change direction 
between the vicinity of the equatorial plane, where they are part of loop 4 and higher latitudes, where they 
are part of loop 3. At t = 132 we found that these radial currents are comparable in magnitude. 

Using the magnitudes of the westward and eastward currents aty = 0, we estimated the magnitudes of 
the different current loops in Figure 5 at the time t = 132, just after the buildup of the SCW from the first 
flow burst. We found that they are similar in magnitude, with 0.36, 0.36, 0.1 6, 0.20, and 0.31 MA for loops 
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1-5, respectively. They are not wirelike but rather 
distributed in sheets or even filling larger space vol- 
umes. We stress that, in estimating the currents of 
loops 1 and 2 (or R1 and R2) from the change in 
cross-tail currents, it is important to include the locally 
closed currents. Using the total eastward and west- 
ward current perturbations in the tail may lead to an 
overestimate of region 1 and 2 currents. 

We also investigated the current closure at the inner 
boundary, taken as a crude proxy of an idealized iono- 
sphere. We demonstrated the closure of the region 
1-type wedge currents through westward currents, as 
postulated in the original SCW picture [ McPherron et 
al. , 1 973], as well as the closure between the region 
1- and region 2-type currents by North-South cur- 
rents. Surprisingly, in a single wedge, the N-S currents 
were larger than the westward currents by a factor 
of almost 2. However, this result should be applied 
with caution to the ionosphere because of our crude 
approximation as a solid plane slab, neglecting finite 
resistance and Hall effects. In particular, the Cowl- 
ing effect [e.g., Cowling, 1932; Baumjohann, 1983; 
Yoshikawa et al., 201 3] should contribute to an enhanced closure through EW currents. 

We should add cautionary remarks on the representation of plasma dynamics by current loops. First of all, 
as demonstrated here, the representation is not unique. Current loops may be combined in different ways. 

If current lines are followed in a simulation (in real space that is not possible), they typically do not close in 
themselves, so investigating current closure may lead one astray, and it depends strongly on whether total 
currents or current perturbations are followed. Finally, none of the representations necessarily implies that 
changing the current in one part of a loop causes a change in the rest of the loop. In a plasma a local change 
in the current may just lead to local closure currents without affecting the global system [Birn etal., 2000]. 

The association of current loops with particular magnetic field changes provides a better picture of cause 
and effect and of the interconnection between the systems. A major role in the substorm dynamics is played 
by the dipolarization in the tail, which includes both an increase of B z in a local equatorial region as well as 
a decrease of \B X \ away from the equator, related to current loops 4 and 1 (and 3) constituting the major 
(region 1 -type) part of the SCW. The additional region 2-type currents are related to two effects, the buildup 
by the azimuthal flow diversion causing increases in \B y \, which also contributes to region 1-type currents, 
and a maintenance by the local increase of pressure (associated with a decrease of B z and current loop 5), 
earthward of the dipolarized region. All of these effects are consequences of the plasma and field dynamics, 
particularly through plasma flows, their braking and diversion, and hence strongly coupled. 

The association of currents with magnetic field changes also brings the modified SCW model closer to 
observations, which are based on magnetic fields rather than currents. The necessity of adding a R2 current 
loop was demonstrated recently by Sergeev et al. [201 1 , 201 4], using Cluster and THEMIS observations in 
connection with geosynchronous observations. The radial width of the westward current associated with 
this loop of ~1 R e also agrees well the width inferred from observations [Sergeev et al., 2014]. However, our 
results differ from the two-loop R1/R2 model of Sergeev etal. [Sergeev et al., 201 1, 2014] by two aspects. The 
westward current in the localized region II (Figure 9) is not necessarily equal to the R2 field-aligned currents, 
due to the presence of loops 4 and 5, which lead to an increase of the westward currents, and loop 2, which 
causes a reduction by radial currents connected to R1 currents. Similarly, the eastward current perturba- 
tions associated with R1 currents are enhanced by the local loop 3 and 4 currents but reduced by loop 2. The 
effects of loop 3 disappear if the cross-tail currents are evaluated from lobe magnetic fields but the effects 
of radial connection with the R2 loop remain. Hence, the R1 and R2 loops inferred by Sergeev et al. [201 1 , 
2014] are not necessarily distinct but may be connected through radial currents and modified through local 
current loops. 


z 



Figure 12 . Schematic of current systems contributing to 
the SCW. 
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